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A modified RSJ-model for the coupling of intrinsic Josephson oscillations and c-axis phonons in the 
high-Tc superconductors Tl2Ba2Ca2Cu30io+5 and Bi2Sr2CaCu20s+5 is developed. This provides 
a very good explanation for recently reported subgap structures in the J-F-characteristic of the c- 
axis transport. It turns out that the voltages of these structures coincide with the eigenfrequencies 
of longitudinal optical phonons, providing a new measurement technique for this quantity. The 
significantly enhanced microwave emission at the subgap structures in both the GHz and THz 
region is discussed. 



The intrinsic Josephson effect between the Cu02 
(multi)-layers in highly anisotropic cuprate supercon- 
ductors has been demonstrated in a characteristic se- 
ries of branches in the current-voltage characteristic in 
c-direction Q. Each branch corresponds to a well de- 
fined number of Josephson junctions in the resistive 
state. Recently the observation of subgap structures in 
the i-U-characteristicof Tl 2 Ba 2 Ca 2 Cu 3 Oio +< 5 (TBCCO) 
and Bi 2 Sr 2 CaCu 2 8+( 5 (BSCCO) has been reported [§,§. 
The structures in the different branches can be traced 
back to the J-V-characteristic of one single Josephson 
junction in the resistive state. These structures seem to 
be an intrinsic effect, as they have been observed in junc- 
tions of different geometries. Furthermore the character- 
istic voltages of the subgap structures are completely in- 
dependent of temperature and external magnetic fields, 
which rules out any relation to the superconducting gap, 
vortex flow or the thermal excitation of quasiparticles. 

In this paper we will show that these subgap struc- 
tures can very well be explained by the coupling of c-axis 
phonons to the Josephson field oscillations in a single re- 
sistive junction. Thereby it turns out that the maxima in 
the I — U-characteristic correspond to the eigenfrequen- 
cies of longitudinal lattice vibrations. A more detailed 
discussion can be found in Q. 

In the RSJ model the total current density through a 
junction 



I/F — j c sin 7 + j qp (E) + D 



(1) 



is the sum of the Josephson current, the quasiparticle 
current and the displacement current density. The (gauge 
invariant) phase difference 7 is related to the electric field 
E in the barrier of thickness b by the Josephson relation 



hj = 2eEb 



(2) 



In the simplest case the quasiparticle current density is 
described by an ohmic law j qp — aE. Here we will use 
a more general form j qp (E) — j qp (E) + cr(E(t) — E) 
with a = j' q p(E) by expanding the field E(t) around the 
time averaged field E. The displacement current density 



D contains the polarization P of the barrier medium, 
D = 6qE + P = eqeE. In the case of high frequency 
Josephson oscillations in the range of phonon frequencies 
it is important to keep the full frequency dependence of 
e(w) or to treat the polarization P as an additional dy- 
namical variable. 

Here we assume that the polarization P is due to a 
c-axis displacement z of ions with charge q, mass M and 
density n in the non superconducting barrier between the 
copper oxide (multi)-layers. For the motion of the ions 
we assume a simple oscillator 



z + u a z + pz = —E, 



(3) 



which is driven by the electric field E in the barrier. In 
this model the contribution of the oscillator to the dielec- 
tric function is given by 



e(uj) = £00 + 



S<4 



lu'q — lo 2 — ipuj 



(4) 



where S = nq 2 /(coMluq) is the oscillator strength of the 
phonon and includes the other contributions to the 
polarizability of the barrier. In this form it is straight- 
forward to generalize the calculation to several phonon 
branches. 

The system is characterized by the transversal op- 
tical phonon frequency ujq, the Josephson plasma fre- 
qency w 2 = 2ebj c / (?ieo e oo) and a characteristic frequency 
u) c = (2e/h)V c , V c being the voltage at the critical cur- 
rent j c . In our case ujj < u>o < u> c and the McCumbcr 
parameter C = (j^/wj is large. 

Using equ. (0) in (Q) and (@) we obtain a coupled set of 
differential equations for the phase 7 and the polarization 
amplitude z, which can be solved numerically. Then the 
time average of the phase velocity (7) = (2e/K)V deter- 
mines the dc voltage V shown in the I —V characteristic. 

The most important features of the model can also 
be understood analytically: Due to the large McCumbcr 
parameter (3 C both the phase and the polarization 



7(*) 
z{t) 



7o 



uit + 71 sin tut, 
z\ cos(ujt + ip) 



(5) 
(6) 
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oscillate primarily with the basic frequency to = (2e/h)V. 
It is also well known that in this case 71 <C 1 and sin (7(4)) 
can be expanded in orders of 71 keeping only the leading 
Fourier components: 

sin 7(f) w sin7o cosut + cos 70 sinwt — —71 sin 70. (7) 

Since the equation of motion for the phonon is linear, 
we can relate the phonon amplitude to the driving force 
7 by the phonon response function of equ. ^ at the 
frequency to. Thus we obtain for the coefficents of sinwt, 
coswt and the dc current i = I/I c : 



= cos 70 + — 0-6171 



or / a 
= sm7 + — 2 £2 H I 7i 



e uj 



1 



i = hp( v ) - ~7i sin 7o , 



(8) 
(9) 
(10) 



where t\, £2 are the real and imaginary part of the phonon 
dielectric function. From these three equations we finally 
obtain for the normalized dc-current 



i{ y)=i qp (V) + -^ 



1 L0 2 , £2 



<~2 



2w 2 V?(w) 



(11) 
(12) 



with e = ei + i(e2 + c/ (weo)) and w = (2e/?i)F. 

This result clearly shows that maxima in the I — V 
characteristic occur near the zeros of the phonon dielec- 
tric function t\ exactly at the eigenfrequencies of longi- 
tudinal optical phonons with wavevector q = 0. Here 
the polarization and the Josephson field oscillations are 
compensating each other, and the energy transfer to the 
phonon is maximal. In contrast to this, the Josephson os- 
cillation amplitude 71 has a minimum at the transversal 
optical eigenfrequency ojq. 

In figures (Q) and (0) this theory is compared with 
experimental results by choosing appropriate values 
for the McCumber parameter (3 C , the phonon frequencies 
a->o,i; their oscillator strength Si and damping pi, and the 
quasi-particle characteristic. We see that the essential 
features are very well reproduced. We also find an ex- 
cellent agreement between the numerical integration and 
the analytical results. Note that both in experiment and 
theory hysteretic regions appear near the phonon reso- 
nances. 

In TBCCO the characteristic voltage is given by V c — 
27.1 mV and the two maxima correspond to LO-phonon 
frequencies ^lo,i = 3.65 THz and i^lo,2 = 4.70 THz. 
The corresponding values for BSCCO are: V c = 21.8 mV, 
Ulo,i = 2.96 THz, ^ LO , 2 = 3.90 THz and i/ LO ,3 = 
5.71 THz. These values are compatible with phonon res- 
onances observed in the infrared reflectivity at 3.7 THz 



and 4.3 THz of polycristalline samples of TBCCO [§. 
Unfortunately for both TBCCO and BSCCO no reliable 
infrared experiments for c-axis phonons in single crys- 
tals are available pQ. Lattice dynamical calculations 
M show phonon modes with large oscillator strength, 
where TIO layers are oscillating against CuO layers at 
somewhat higher frequencies. The dispersion of these 
LO-phonons in c-direction is small, which causes a small, 
but finite damping p in the model described above. 
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FIG. 1. Experimental, analytical and numerical J-F-cha- 
racteristic of TBCCO for f3 c = 375, w ,i = 3.28 thz, 
Si = 2.87, pi = 0.2 THz, u ,2 = 4.45 THz, & = 0.68, 
p 2 = 0.4 THz, £oo = 8.45. 
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FIG. 2. Experimental, analytical and numerical J-F-cha- 
racteristic of BSCCO for j3 c = 800, w ,i = 2.80 THz, 
Si = 1.35, pi = 0.08 THz, w ,2 = 3.62 THz, S 2 = 1.5, 
p 2 = 0.05 THz, ujo,3 = 4.91 THz, S3 = 1.65, ps = 0.2 THz, 
too = 7.5. 

A similar calculation can be used to determine the mi- 
crowave power radiated by Josephson oscillations: the 
emission is enhanced by about one order of magnitude 
near the maxima of the subgap structures, suggesting 
these voltages for an experimental attempt to demon- 
strate the microwave emission of the intrinsic Josephson 
oscillations directly. Also the observed enhancement of 
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the microwave emission in the GHz region, if the junction 
is biased near the subgap structures, can be explained: 
the vanishing differential conductivity near the maxima 
of the I — V-cmve leads to an amplification of voltage 
fluctuations caused by thermal fluctuations of the cur- 
rent. 

Further theoretical investigations suggest that also the 
mixing of two external frequencies (in the THz region) 
down to the difference frequency (GHz) is much more 
efficient, if the external frequencies are near the maxima 
of the subgap structures. 

The modifications of the intrinsic Josephson oscilla- 
tions through the dielectric constant of the phonons is in 
complete formal analogy with the coupling of a macro- 
scopic Josephson junction to an external electromagnetic 
system with a resonance, i.e. a cavity. In such systems 
hysteretic structures in the /-^-characteristic, which cor- 
respond to the subgap structures considered in this pub- 
lication, have been reported repeatedly (|]. 

To summarize, the recently discovered subgap struc- 
tures in the /-^-characteristic in the intrinsic Josephson- 
effect in high-T c -superconductors are explained by a cou- 
pling of Josephson oscillations to optical c-axis phonons 
within a modified RSJ-model. Apart from perfect repro- 
duction of the experimental data for both BSCCO and 
TBCCO new insights in the physical interpretation of 
the peak structures are obtained. The peak position is 
identified with the eigenfrequency of the longitudinal op- 
tical phonon, which provides a new possibility of a direct 
measurement of this quantity, which is normally hard to 
determine in optical experiments. Significantly enhanced 
microwave emission near the maxima of the subgap struc- 
ture in the GHz region is explained and theoretically pre- 
dicted in the THz region. 
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